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1. Introduction 

Domain-wall simulations use lattices with a finite number of points N s in the 5th dimension 
[j|], and so a breaking of chiral symmetry occurs. Only in the theoretical limit in which N s = °° the 
chiral modes can fully decouple from each other, yielding an exact chiral symmetry Here we study 
these chiral violations using perturbative calculations and computing three quantities: the residual 
mass m res , the difference A = Zy — Z4, and c m \ x , a chirally-forbidden mixing (which is then nonzero 
at finite N s ) of an operator which measures the lowest moment of the g2 structure function. 

We have studied the dependence of these three quantities on N s and the domain-wall height 
M, and calculated the deviations from the N s = °° results when N s is limited to small values, of 
O(10). We have hence repeated the computations for several choices of N s and M. A thorough 
exploration of large regions in the two-dimensional space spanned by N s and M would be instead 
quite expensive for Monte Carlo simulations, and perturbation theory seems the more practical and 
cheaper way to gather hints of what is happening when the parameters are moved in this space. 

In order to carry out these calculations one must use the Feynman rules which correspond to 
the theory truncated at finite N s , and thus we also had to compute the required propagator functions. 

We have calculated the same quantities with the plaquette action ^ as well as with improved 
gauge actions pi], since in numerical simulations it was observed that these improved gauge actions 
(especially DBW2) reduce the chiral violations. We refer to [Q, [|] for the actions, notations and 
conventions used, and in particular for the expressions of the domain-wall fermion propagators at 
finite N s . Here we only remind that < M < 2 and that the chiral projectors are P± = (1 ± 7§)/2. 

This domain- wall formulation [Q] corresponds to having several flavors of lattice Dirac fermions 
which are mixed via a mass matrix in a very special way, so that a large mass hierarchy is gener- 
ated. To determine the chiral modes one must diagonalize (in the fifth dimension) this mass matrix, 
which however is not hermitian. Its square must then be considered, which means the second-order 
operators DD^ and D^D. They are hermitian and nonnegative and give a well-behaved spectrum. 

A rotation of the 5-dimensional quark fields \\f s {x) to the basis which diagonalizes the mass 
matrix gives finally the expression of the chiral mode: 

Xo(x) = £(/»+<->*(*) +P-w£'->,(x)), 

where from now on we put wq = 1 — M. We can see from the damping factors Wq _1 and Wq'~ s that 
the chiral mode is exponentially localized near the two walls at s = 1 and s = N s . However, the 
domain-wall height M, which is not protected by chiral symmetry, undergoes an additive renormal- 
ization, so that wo is also additively renormalized. 

The standard chiral mode used in Monte Carlo simulations is then not Xo( x )- ^ contains instead 
only the quark fields exactly located at the boundaries: 

q{x) = P+ Yi (x) + P yn s 0) , q(x) = Wn s i x ) p + + Wi i x ) p - ■ 

These physical quark fields q(x) are more convenient to use than %q{x): they do not contain wq and 
avoid the problem of its renormalization. 

At finite N s an additional issue arises: Xo( x ) itself is no longer the exact expression of the 
chiral mode. In fact, Xo{ x ) at finite N s is an eigenvector of the mass matrix only up to terms of 
order N s e~ Ns<x ( ', where ce(0) is a constant determined by 2cosh(a(0)) = (1 + Wq)/|wo|- 
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2. Residual mass at tree level 

The calculation of the propagator of the (approximate) chiral fields q(x) gives (for m = 0) 
% sin p M (1 - e - 2N * a ^) + e - N ° a ^ ■ 2W(p) sinh(a(/?)) 



(?(-p)?(p)> 



1 -W(p)e a W -^"^(l -W(p)r a W) 



where W(p) = 1 -M + 2£ A sin 2 ^ and 2cosh(a(p)) = (1 +W 2 (p) + Ia sin 2 /? A )/|W(/?)|. In 
the limit of small momentum this 4-dimensional propagator becomes 



# + Hff'(l-wg) 

p«i " p 2 + wf s (\-wf 2 ' 



(q(-p)qip)) =-(1-^) 



We can thus see that, although in the bare Lagrangian all quark fields are massless, the truncation 
of domain-wall fermions at finite N s generates already at the tree level a nonvanishing residual 
mass of the physical fields: am^] = — w^ ! (l — Wq) = —(1 —M) Ns M(2 — M). As expected, this 
tree-level residual mass vanishes when N s becomes infinite 1 . Its sign can be inferred from the 
general expression of a fermion propagator of mass jl for small momentum in Euclidean space: 
(— if 1 + (p 2 + pL 2 ) = \ / {ip 1 + }l) . Since we work with even N s (where the fermion determinant 
can be proven to be positive), m^ s is always a negative quantity. With our calculations we have 
thus reproduced, up to a sign, the result for mPe S found in |5[ ^, [7| H], where it was derived by 
considering the quadratic operator D^D, which could perhaps explain the sign discrepancy. 

3. Physical propagator at one loop 

At one loop we can write 2 



(q{-p)q{p))l loop = J- + ~ if- T MP) -., a, , . 

ip-w iP-Wq\\-w z q) ip-w Q \\-w z Q ) 

1-w 2 



^-^(l-w 2 )-(l-w 2 )I,(p)- 
The general form of E q (p) for m = is (where we call for brevity g 2 = (g 2 ) /\6n 2 )Cp) 



Up) = -+V 47 ^logaV+Si - (¥-<(i -<)) f^i3 

1 — Wq I a V 1 / v 1— Wq 





'indeed, since wo = e~ a (°\ it is easy to see that the terms which are proportional to Wq° = e~ Ns<x ^ rapidly approach 
zero when N s becomes large. 

2 Evidencing the damping factors of the external legs and the loop integral L st (p), the structure of is 

v I \ r*N s tN s 1 \l N,-s 2N, —(N s —s)\ D / s _l 2N, —(s — 1)\ r, 
MP) = E,=l I f =l " 2AT K - w w ) P + + K - w w ) P - 

l-w 

-W() j IK - w V ) P + + K ~< V ) P - 

l-w v 

r / \ 1 IV W,-t 2iV, -{N,-t)\ D . / t-1 22V. -(f-l)\ D /tin 

2WT K ~ w o w o ) p - + { w o ~ w o V ) p + ( 3J ) 

l-w 

-w {{W -w„ w )P- + {w -w w )/Vj j— ^5 J- 
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The most important difference in L q from its expression at infinite N s is the appearance of a totally 
new contribution, Zo> proportional to I /a and associated with the breaking of chiral symmetry. Zo 
comes from the terms of l, st (p) which are of order zero in p, and acts as a mass correction term. 
Since 

1 — 1 — Wq 

m-pmp)h loop = - — — J—- — = — — — ^yZ w , 

the one-loop radiatively induced mass is given by 

««« = -w ! (l - W ) -g Z . 

Thus, £() generates a finite additive renormalization to the residual mass when N s is not infinite 3 . 
The factor Z w = 1 — 2g 2 zZ^ Wq/(1 — Wq) = l+g 2 z w generates the additive renormalization to w>o at 
this order as we can see from (1 — wfyZ w = 1 — (wq + g 2 Lj,) 2 + 0(g A ). 

The renormalization of a composite operator q{x) Oq{x) which is multiplicatively renormaliz- 
able can also be expressed in a simple way: 

{{qOq)qq)\ loop = - jj- T • A (p) • ^- -, 

where Ao(p) contains also the contribution of the damping factors, and takes the form 

Ao(p)=S 2 (-Y { o s ' M hogaV+Bo 

for a logarithmically divergent operator. The coefficients of the divergences turn out to be different 
from their continuum values, and they depend on N s and M. It is only when N s = °° that the 
anomalous dimensions become the ones calculated in the continuum theory. In particular, the 
vector and axial- vector currents acquire a nonzero anomalous dimension at any finite N s : 

(N sM ) _ lw 2N s (, N 2N S l-H'Q \ ( N U 1 \ 7 W l \ 

Furthermore, the residual mass as well as renormalization factors and mixing coefficients turn out 
to lose gauge invariance when N s is not infinite. Although numerically the deviations from gauge 
invariance remain in most practical cases rather small, this is another of the pathological features 
of the domain- wall theory truncated at finite N s . It could be that this is a limitation of perturbation 
theory, but it could be that a small gauge dependence is also present in numerical simulations. 

These pathologies could actually be related to the mismatch between the (simplified) chiral 
modes which are actually used, and the true chiral modes (the ones which contain wq). Notice that 
at finite N s there is an additional mismatch, because terms of order Af s e~ NsCt ^ and higher, which are 
present in the true chiral modes for < °°, are here missing as well. Thus, if calculations with the 
true chiral modes would be gauge invariant and reproduce the continuum anomalous dimensions, 
the missing pieces from these mismatches could then account for the above pathologies. 

At one loop two diagrams contribute to Lq and so enter in the calculation of the residual mass: 
the half-circle (or sunset) and the tadpole diagrams. We have automated the calculations of the 
half-circle diagram (as well as the vertex diagrams for A and c m ; x ) by developing suitable FORM 
codes [|l0|], integrating afterwards the corresponding expressions by means of Fortran codes. With 
these programs we are able to compute matrix elements for general values of Af s and M. 

3 Of course higher loops and nonperturbative effects give further contributions to the shift of the residual mass. 



4 



Chiral violations from one-loop domain wall fermions 



4. The tadpoles 

The behavior of the tadpole diagrams as N s and M change is particularly important. 

The tadpoles do not contain pure 5-dimensional quark propagators, and so for them the L st (p) 
of Eq. ( jO| ) is diagonal in the fifth-dimensional index, and also proportional to ($—Ar/a) G^ v (k), 
where G^ v is the gluon propagator. This is the same integrand of the tadpoles for Wilson fermions, 
where in the case of the tadpole diagram contributing to Ei it gives the result (in a general covariant 



gauge) 7} = 87T 2 Z (1-1/4(1- A)), with Z = 0.154933390231 ... a well-known integral Q. 

It is then clear that for domain-wall fermions the behavior of the tadpole diagrams as a function 
of N s (and M) is completely determined by the damping factors in the fifth dimension, (see Eq. 
(j3~l|)). Their general effect can be already seen by looking at their leading contributions for large 
N s . In this approximation the damping factors enter the game in either of these combinations: 



(wo) 



Es-1 N.—s at N.—l V / 2\s-l V / 2\N s -s 

.v=l s=l s=l 



These are indeed the leading expressions, in units of 7^ = (1 — Wq)Ti /(l — w Ns ) 2 , for the tadpole 
contributions to (respectively) £o and £i, in the limit of large N s . Already from these asymptotic 
expressions (before computing the exact results) we can immediately see that the tadpole of £0 
vanishes when N s = 00, while the tadpole of £1 gives in this limit the known Wilson number, 7}. 

Thus, the damping factors play a primary role in determining the values of the domain-wall 
tadpoles. After calculating their exact expressions, which include all subleading terms in N s , the 
tadpole contribution to £ turns out to be equal to 

1 2N S 

4T d N s (l+w y ')w s -2w Q s — 

l-w . 

while the tadpole contribution to £1 turns out to be equal to 



1 2N S 

(l+w^ +1) )^V- 2 ^wf 
1 — w. 



"o 

The values of these tadpoles present wide variations with N s and M, so that sometimes they turn 
out to be small while in other situations they become large. This suggests that some care should be 
used when talking about tadpole dominance in relation to domain-wall fermions. It also happens 
that both tadpoles (of £ an d £1) even decrease toward zero for M — ► or M — > 2. 
A central point is that there are two kinds of tadpoles in the game here: 

• the tadpole of order zero in p, which tends to zero for N s — * °°, and which contributes to £0 
and the residual mass; 

• the tadpole of order ap, which tends to its Wilson value for N s — > °°, and which contributes 
to £1 and the renormalization factors. 

They behave quite differently, and this is because the ip 1 of the first order of the self-energy flips 
the chirality of some damping factors, which then combine in a different way. We stress that this is 
quite unlike the Wilson case, where the tadpole of £0 is just proportional to the tadpole of £1 : 

rp(Wilson) . rp(Wilson) 

7 (Z ) -~ 47 (£i) • 
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Tadpole improvement seems then not to be appropriate for the residual mass: the tadpole 
which contributes to m res goes to zero for large N s or for M — > 1, and for small N s it assumes a 
wide spectrum of values. Our interest is small N s , where it is unclear what tadpole improvement 
(or resummations) could mean. Moreover, no tadpole enters at all in the calculations of A and c m &. 

For large N s the tadpole of Ei is rather close to its Wilson value, and that is why tadpole 
improvement could be used in the calculations of the renormalization factors (in the large N s limit). 

5. Residual mass at one loop 

Our one-loop perturbative calculations show that the numerical deviations from the case of 
infinite N s depend, apart from N s (and to a smaller extent from go), very strongly on the choice of 
M. We can observe that the deviations from the case of exact chiral symmetry are rather pronounced 
when M ~ 0.1 or M ~ 1.9. The values of arrives turn out to be positive only for M > 1.2 (at least 
for even N s and if the coupling is not very small), and our results suggest that the minimal amount 
of chiral violations is attained for M ~ 1.2. This is then the optimal choice of M from the point of 
view of one-loop calculations, and corresponds to the renormalization of M, which is not protected 
by chiral symmetry and is then moved by radiative corrections away from its free field value M = 1 . 
One can conjecture that higher-loop corrections and nonperturbative effects would shift this optimal 
value further on, until the minimal point is eventually reached around M ~ 1.8 (which provides the 
smallest residual mass in Monte Carlo simulations). 

We also can observe that for M = 1.9 the residual mass at N s = 12 is larger than at N s = 8, 
and at A^ = 16 is even larger. For a detailed discussion of these phenomena which occur near the 
borders of the allowed values for M we refer to [g]. 

With improved gauge actions we can still see that the residual mass arrives is positive only for 
M > 1.2. This also shows that improved gauge actions do not behave too differently in terms of 
the additive renormalization undergone by wo- Employing improved gauge actions produces, not 
surprisingly, a suppression of m res when one carries out the comparisons at the same value of the 
coupling. The Iwasaki action gives a stronger suppression than the Liischer-Weisz action, and there 
seems to be a monotonic decrease of the residual mass as c\ grows. The DBW2 action is indeed 
the most effective in generating large suppressions. 

If comparisons between the various actions are instead made at the same energy scale, the 
picture that comes out is different from naive expectations. For example, for quenched QCD at 2 
GeV one has to take j3 = 5.7 for the Liischer-Weisz action, j8 = 2.6 for the Iwasaki action, and 
j8 = 1.04 for the DBW2 action 4 . The m res numbers for the Iwasaki action are then rather close to 
those of the DBW2 action, and surprisingly they lie in general slightly above the plaquette values. 

However, for the quenched DBW2 action at 2 GeV one has = 5.77, which is rather large, 
and so the one loop results cannot perhaps be trusted so easily 5 . Moreover, these values of the 
couplings are determined from numerical simulations, and they then contain informations of a 
nonperturbative nature, so that a mismatch can arise when one only takes into account the results 
of the one-loop diagrams calculated for these values of the couplings. 

Many numerical results, which we cannot include here for lack of space, can be found in B]|. 

4 We use /3 = 6/gg also for improved actions, instead of /3' = 6(1 — 8c'i )/g\- 

5 In this case the Liischer-Weisz action gives the largest m res suppression, and indeed at 2 GeV is still close to 1. 
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6. Bilinear differences, and a power-divergent mixing 

Since Zy ^ when chiral symmetry is broken, the difference between these renormalization 
constants, A = Zy — = — (Z$ — Zp)/2, provides an estimate of chirality-breaking effects. 

The amount of chirality breaking connected to A follows a pattern similar to the one of the 
residual mass: A is rather large for small N s or |1 — M\ ~ 1, it decreases when N s grows or when 
|1 — M\ tends towards zero, and the violations of gauge invariance are very small. 

The numbers for A come out much smaller, at a given M and N s , than the ones for the residual 



mass. For quantities such as four-quark operators it was suggested in [|12|, |13|] that their chiral 
violations are of 0(nr} es ). Given the smallness of the numbers that we have obtained for A, it is 
possible that something similar is also occurring here. 

We have also calculated the mixing of the antisymmetric operator 

Od, =q{x)y [A y 5 D l] q{x) 

with an operator of lower dimension, 

c mix - l -q{x)o A \y 5 q(x). 

The operator Oj, enters in the calculation of the first moment of the g2 structure function, and has 
been simulated using quenched domain-wall fermions with the DB W2 gauge action [[14]] . 

The power-divergent mixing of O l j l on the lattice is only caused by the breaking of chirality, 
and hence it provides a quantitative measure of chiral violations. In the theoretical limit N s = °° 
one has instead c m & = and so 0^, becomes multiplicatively renormalized. 

The chiral violations associated with c m « are rather small, and thus they also seem to be of 
higher order in m res . The pattern of the deviations from the case of exact chirality is the usual one. 
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